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Abstract 

We study the approximation properties of the class of nonstation¬ 
ary refinable ripplets introduced in m- These functions are solution 
of an infinite set of nonstationary refinable equations and are defined 
through sequences of scaling masks that have an explicit expression. 
Moreover, they are variation-diminishing and highly localized in the 
scale-time plane, properties that make them particularly attractive 
in applications. Here, we prove that they enjoy Strang-Fix condi¬ 
tions and convolution and differentiation rules and that they are bell¬ 
shaped. Then, we construct the corresponding minimally supported 
nonstationary prewavelets and give an iterative algorithm to evaluate 
the prewavelet masks. Finally, we give a procedure to construct the 
associated nonstationary biorthogonal bases and filters to be used in 
efficient decomposition and reconstruction algorithms. 

As an example, we calculate the prewavelet masks and the non¬ 
stationary biorthogonal filter pairs corresponding to the nonsta¬ 
tionary scaling functions in the class and construct the corresponding 
prewavelets and biorthogonal bases. A simple test showing their good 
performances in the analysis of a spike-like signal is also presented. 
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1 Introduction 


A ripplet is a function / whose integer translate are totally positive [26] . 
i.e. for any ordered real numbers xi < ■ ■ ■ < Xr, and any ordered integers 
ai < ■ ■ ■ < Uri r > 1/it holds 

det {f{xi - ai)) > 0. (1.1) 

Total positivity implies that the integer translates of / are variation dimin¬ 
ishing, i.e. for any hnite sequence c = {ca} 

S-(5^c„/(--a)) <S-(c), (1.2) 

where S~ denotes the strict sign changes of its argument. The disequality 
fll.21) in turn implies that the system {/(• — a)} has shape-preserving prop¬ 
erties, which are known to play a crucial role in several applications, from 
approximation of data to CAGD [TT[ |28] . 

The concept of a ripplet was hrst introduced by Goodman and Micchelli 
in 

[T5] . where the authors focused their interest on two-scale refinable ripplets, 
i.e. ripplets that are solution of a two-scale refinable equation 

ip = ^ -a), (1.3) 

OL 

where the scaling mask a = {oq} is a suitable real sequence. Well known 
examples of rehnable ripplets are the cardinal B-splines, i.e. the polynomial 
B-splines on integer nodes. Starting from the seminal paper of Goodman 
and Micchelli, many families of two-scale rehnable ripplets were constructed 
(see, for instance, 0. ini)- More recently, M-scale rehnable ripplets, with 
dilation M greater than 2, were addressed in [16] and rehnable ripplets with 
dilation 3 were constructed in [20] (see also [2T]h 

The interest in rehnable ripplets lies in the fact that they give rise to con¬ 
vergent subdivision algorithms for the reconstruction of curves and the limit 
curves they generate preserve the shape of the initial data [T3l |28] . Rehnable 
ripplets can also be proved to solve the cardinal interpolation problem [28] . 
For instance, the construction of cardinal interpolants by by means of the 
rehnable ripplets in d was addressed in [30] . 
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Refinable ripplets have good properties not only in the context of geomet¬ 
ric modeling and function approximation but they also enjoy some optimal¬ 
ity properties useful in signal processing. In fact, rehnable ripplets induces a 
multiresolution analysis in L 2 ( 1 K) that allows us to generate a nested sequence 
of wavelet spaces. Actually, it is always possible to construct compactly sup¬ 
ported semiorthogonal wavelet bases starting from a rehnable ripplet m- 
Moreover, rehnable ripplets have asymptotically, i.e. when their smoothness 
tends to inhnity, the same optimal time-frequency window achieved by the 
Gaussian function [2]. Since the rate of convergence can be proved to be 
very fast for a large class of ripplets, including, for instance, the rehnable 
ripplets introduced in na, rehnable ripplets can approximate the Gaussian 
with high accuracy giving rise to efficient algorithms for signal analysis [T]. 
Finally, a ripplet can be seen as a discrete kernel satisfying a causality prop¬ 
erty so making the rehnable ripplets particularly attractive in the scale-time 
analysis of signals [T2] . 

All the rehnable ripplets quoted above are stationary in the sense that 
they satisfy the functional equation fll.dp with the same mask sequence at 
each dyadic scale. For this reason usually fll.dp is referred to as a stationary 
two-scale equation. From the point of view of signal processing this means 
that the same analysis and synthesis hlters are used at all dyadic scales [3T]. 
Nevertheless, the use of the same set of hlters at each scale does not give 
great hexibility in applications, especially when some preprocessing steps 
with diherent hlters are required. From the functional point of view, the 
use of diherent hlters at diherent scales gives rise to a nonstationary mul¬ 
tiresolution analysis (see, for instance, i, la, 1221. 1231). A nonstationary 
multiresolution analysis can be generated by a set of nonstationary refinable 
functions, i.e. an inhnite set of functions { 93 ™ : m G Z+} satisfying an inhnite 
set of nonstationary two-scale equations 

= ^ aff m e Z+ , (1.4) 

for some inhnite sequence of masks {a™ : m G Z+j, each mask a™ = {ajj*} 
being diherent at each dyadic scale. 

The use of diherent scaling masks at diherent scales allows us to con¬ 
struct rehnable functions endowed with properties that cannot be achieved 
in the stationary setting. For instance, the nonstationary process generated 
by the B-spline masks of increasing support gives rise to a rehnable func¬ 
tion that is compactly supported while belongs to G°°(M) [IHl EH]- Other 
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families of nonstationary refinable functions based on pseudo-spline 

masks were constructed in [21j. Exponential splines too can be associated 
to a nonstationary process ii- They reproduce exponential polynomials 
and the associated wavelet bases can be successfully used in the analysis of 
signals with exponential behavior |32]. Families of nonstationary rehnable 
ripplets were introduced in [S] and [12]. In particular, the latter are highly 
localized in the scale-time plane, a property that is crucial in many appli¬ 
cations. For this reason, in the present paper we focus our interest in this 
family of nonstationary rehnable ripplets and prove that these ripplets en¬ 
joy several properties which are relevant in the context of both geometric 
modeling and signal processing. Then, we construct the associated wavelet 
bases. We notice that, since the rehnable ripplets we are considering are non 
orthogonal, compactly supported orthogonal wavelets belonging to the space 
generated by their translates do not exists. This motivates us to construct 
nonstationary semiorthogonal prewavelets. In fact, giving up the orthogonal¬ 
ity condition we can build wavelet bases with compact support. Moreover, 
we give a procedure to construct the nonstationary compactly supported 
biorthogonal bases associated with the nonstationary rehnable ripplets we 
are considering. Since we are interested in implementing efficient nonsta¬ 
tionary decomposition and reconstruction algorithms for signal processing, 
we construct also the corresponding pairs of nonstationary hlters to be used 
in the analysis and synthesis of a given discrete signal. 

The paper is organized as follows. In Section 2 we give some basic deh- 
nition concerning nonstationary multiresolution analysis and wavelet spaces 
and recall some results about the existence of solutions of the nonstationary 
rehnable equations fll.4p . The class of nonstationary rehnable functions we 
are interested in is described in Section 3, while in Section 4 we proved some 
approximation properties that were not addressed in [12]. In Section 5 we 
construct the nonstationary prewavelet bases associated with the nonstation¬ 
ary rehnable ripplets in the class and give an efficient algorithm to evaluate 
nonstationary prewavelet masks. The construction of compactly supported 
biorthogonal bases and hlters, which give rise to efficient decomposition and 
reconstruction formulas for discrete signals, is addressed in Section 6. Finally, 
in Section 7 some examples of nonstationary masks and rehnable bases are 
given and the corresponding nonstationary hlters are constructed. A simple 
test on the analysis of a spike-like signal is also shown. 
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2 Nonstationary Multiresolution Analysis and 
Wavelet Spaces 

Wavelet spaces are usually constructed starting from a multiresolution anal¬ 
ysis that is a sequence {V^} of nested subspaces which are dense in L 2 (M) 
and enjoy the separation property. Thus, the corresponding wavelet space 
hh™ C is dehned as the orthogonal complement of Id™ in [7]. 

In the stationary case all the spaces {V™} are generated by the dilates and 
translates of a unique rehnable function (p E V^, solution of the stationary 
two-scale equation fll.dp . while the spaces W™ are generated by the dilates 
and translates of a unique wavelet tjj G C V^. 

In contrast, in the nonstationary setting any space V™ (resp. IT™) is 
generated by the 2“™-shifts of a different rehnable function (p™ (resp. wavelet 
-0™), which are not dilates of one another. As a consequence, the spaces {T™} 
and {IT™} are not scaled versions of the spaces T° and IT*’, respectively. 

Nonstationary multiresolution analysis are addressed in several papers in 
the literature (see, for instance, i, ffl, i, la, ra, ra, ra and references 
therein). Here, following [8] and [29], we say that a space sequence {T™ : 
m G Z+} forms a nonstationary multiresolution analysis of L 2 (ffi) if 


(z) T™ C T™+\ m G Z+; = i^ 2 (M); (m) = {0} 


{iv) for any m G Z+, there exists a L 2 (M)-stable basis ip^ in T™ such that 
T™ = {(p™(- - 2-™a), a G Z} . 

Property (z) implies that the refinable functions {(p™ : m G Z+} satisfy a 
set of nonstationary refinable equations, i.e. 


^ a™ (p™+^(--2-(™+i)a), mGZ+, (2.1) 


m+1 


for some sequence of scaling masks {a™ : m G Z+}, where a™ = (a™ : a G 
Z} G i 2 {'^) and 



( 2 . 2 ) 


Properties (zz) and (zzz) are always true if any <p™ is compactly supported. 
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As for (iv), the function system — 2“™'a) : a G Z}, m G Z+, is L 2 (ffi)- 

stable if and only if the Fourier transform of has no 2™'+^7r-periodic real 
zeros (cf. [8]). 

The existence of a unique set of functions : m G Z+} solution to 
fl2.ip . as well as their properties, are related to the properties of the mask 
sequence {a™ : m G Z+}. Since we are interested in the the case of compactly 
supported masks, we assume that any mask a™ is compactly supported with 

supp (a™) C 12 c Z, m G Z+ , (2.3) 

and there exists a fundamental mask a = {a^ : a G 12} satisfying the sum 
rules 

02(1+7 ~ 1 > 7 ^ ^) (2-4) 

such that 

\aff — Ual < oo , a G 12, (2.5) 

(cf. [S]). Even if the assumption above excludes the mask sequence associ¬ 
ated with the up function and the nonstationary mask sequence considered in 
[22] . nevertheless it covers many nonstationary scaling masks, such us those 
ones associated with the exponential splines 0. ra and the mask families 
associated with the ripplets introduced in [5], [19]. For the case of mask 
sequences with growing support we refer the reader to [1], [29] . 

The nonstationary rehnable equations fl2.ip can be associated to a non¬ 
stationary cascade algorithm [H], which generates at any iteration fc G Z+ 
the sequence of functions {h™ : m G Z+} by 

'Cl = E < e ™ e ■ (2-6) 

aGO 

Without loss of generality, we assume that the starting function is the 
same for all m G Z+, i.e. h™ = ho, where hg is a given L 2 (M)-stable function 
with ho(0) = 1, so that h^(0) = 1 for any /c,m G Z+. 

The convergence of the cascade algorithm is related to the spectral prop¬ 
erties of the fundamental transition operator T : dehned as 

(T A)a = 2 d 2 a -/3 , a E h, (2.7) 

/3gZ 
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where 

a = {tta e Z} , Cla = 0/3 0/3-0 ) tt G Z, (2.8) 

/3gZ 

is the autocorrelation of the fundamental mask a. The cascade sequence 
{h™ : fc G Z+} converges strongly to in L 2 (M) as k ^ oo, uniformly in m, 
if and only if the fundamental transition operator has unit spectral radius, 1 
is the unique eigenvalue on the unit circle and is simple [HI Th. 1.3]. Under 
these hypotheses the sequence {ip"^ : m G Z+} converges strongly to the 
solution of the stationary rehnable equation 

ip = 2 '^ aaip{2- -a). (2.9) 

We notice that any nonstationary mask sequence {a™ : m G Z+} having 
a B-spline mask as fundamental mask, gives rise to a convergent cascade 
algorithm PH. 

The nonstationary rehnable equations in the Fourier space read 

^(w) = kl”*(e"*^)(^+^(a;), m G Z+, (2.10) 

where the Laurent polynomials 

mGZ+, zeC, (2.11) 

oGO 

are the mask symbols. We notice that any rehnable function ip^ is normalized 
so that ^(0) = 1. In case the nonstationary cascade algorithm converges, 
the Fourier transform of (p'^ is given by 

00 

^{u) = Yl , wgM. (2.12) 

k=m 

Given a nonstationary multiresolution analysis {U™ : m G Z_|_}, at any 
level m G Z+ the wavelet space hF™' is dehned as the orthogonal complement 
of U™' in i.e. 

lU™ = © U'" , m G Z+ . (2.13) 

Any wavelet space is generated by the 2“™'-shifts of a wavelet function '0™, 

i.e. 

= spm {^^{■-2-^a),aeZ} , (2.14) 
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any rjj'^ being different at any scale m. The existence of a set of generating 
wavelets G Z+} is always assnred when all the rehnable fnnctions 

G are compactly snpported. Moreover, it is always possible to 
constrnct compactly snpported L 2 (K)-stable wave-lets associated with com¬ 
pactly snpported L 2 (K)-stable rehnable fnnctions [H §4], 


3 A Class of Bell-shaped Nonstationary Re- 
finable Ripplets 


The class of nnivariate compactly snpported nonstationary masks and rehn¬ 
able fnnctions we are interested in was introdnced in [19]. Let ns denoted 
the masks in the class by 

{a<"'’“>:mez+}, a'”--"'= , a^T’} , m e Z+, (3.1) 


where n is an integer > 2, related to the snpport of the mask 
The explicit expression of is as follows. For any n, the starting mask 

a(">o) Pas entries 




= = 5 


1 

2 ’ 


“0 

_ Q ^ otherwise , 
while the entries of the higher level masks m > 1, have expression 


(3.2) 


(n,m) 

a 


1 


n+l 
a 


+ 4 ( 2 ”^''^ 



0 < a < n -I- 1, 


a 


0 , 


otherwise. 


(3.3) 


(We assnme (”) = 0 when a < 0 or a > n.) 

We notice that /i > 1 is a real parameter that acts as a tension parameter. 
In fact, the larger /i the faster the factor 2”* ^ in the nonstationary mask 
coefficients goes to 1 when m —)■ cx). This means that the nonstationary 
process behaves in practice as a stationary one if we choose large valnes of /i. 
Thns, the more interesting cases are obtained for valnes of fi that are close 
to 1. 

For any n and /r, the mask is compactly snpported on [0, n -|- 1] and 

is hell-shaped^ i.e. its entries are positive, centrally symmetric and strictly 





increasing on [O, 

We notice that for any n >2 the 0-level scaling mask is the mask of the 
characteristic function of the interval [0,1], while the m-level scaling masks 
al"’™'), m > 0, are related to the class of stationary masks introduced in ra- 
Any mask sequence : m E Z+} can be associated with the set of 

nonstationary rehnable equations 


A 


meZ+, (3.4) 


aGa, 


{n,m) 


where 


^(n,m) _ cni-m ^ ^ ’ 

^ ^ ( [0, n -I- 1], m > 0 . 


(T^a = supp 


(3.6) 


From [19] it follows that for any n > 2 and p > 1, m E Z+, is 

compactly supported with 


SUpp = [O, L(^n,m)] 


[o, f + l] , m = 0 , 

[0,2-™(n +1)], m>0. 


(3.6) 


and belongs to C” ^(K). Moreover, any system 


^(n.m) ^ {^(n,n.)(. _ 2-"'a) : a € Z} , 


(3.7) 


is linearly independent and L 2 (M)-stable, forms a partition of unity, is totally 
positive, and enjoys the variation diminishing property. 

For any n >2 and /i > 1, the mask sequence : m E Z+} has the 

sequence 


a(”) = {al">,0<a<n + l}. a'”) = ^. (3.8) 

as fundamental mask [19]. Since is the mask of the B-spline of degree n 
having integer knots on [0, ri-1-1], the sequence of functions : k E Z+}, 

generated by the cascade algorithm 

^ h^”^+')(--2-(™+i)«), meZ+, (3.9) 

(n,m) 
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converges strongly to G L 2 (M). The convergence of the cascade algo¬ 
rithm implies that the Fonrier transform of is 

OO 

^(n,m)(^) _ Yl , ueR, (3.10) 

k=m 

where 

^ . (3.11) 

A straightforward computation gives 

= yi + j) 

(1 + + 2(2’+"''' - 1) .- + 1) ^ (3.12) 

= (;2) = yl(0) yl(rx-l,m) ^ m > 0 . 

It is worthwhile to observe that any symbol (z) is a Hurwitz polynomial 

m, i-e. it has only zeros with negative real part. Moreover, the fundamental 
symbol of A^”'’™'^(z) is the B-spline symbol 

n+l .. 

At’'>(z) = 5^ 4"! = — (1 + ^)"+i. (3.13) 

a=0 


4 Properties of the nonstationary refinable 
ripplets 


In this section we analyze some properties of the rehnable ripplets that 

are useful in both geometric modeling and signal processing applications. 

First of all, let us denote by 5the 2^-dilates of the B-spline of degree 
n with knots on 2“™'Z, normalized so that i?*^"’’”*)(0) = 1, n > 1, m G Z+. In 
the Fourier space 5satishes the rehnable equation 


5(n,m)(^) ^ ^(n) 

Its Fourier transform is given by 

OO 

5(n,m)(^) _ J-J ^(n) 


k=m 







The refinable functions are generated by a convolution low involving 

the B-spline 


Theorem 4.1. Forn > 3, m G satisfies the convolution property 

^ ^(0,1) ^ ^(n-1,0) ^ 

^in,m) ^ ^ ^ m > 0 . 

Proof. From fl3.10p . fl3.12p and fl4.2p it follows 


(4.3) 




J] A 


k=m 


g ^2^+T 




k=m 




for m > 0, and for m = 0. The claim 

follows by applying the inverse Fourier transform to the relations above. □ 


As a hrst consequence of the theorem above, we can prove that satishes 

suitable Strang-Fix conditions. 


Corollary 4.2. For any n > 2 andm > 0, has a zero of order n — 1 for 

oj = 2™'+^7ra, a G Z\{0}. has a simple zero for u = 27ra, a G Z\{0}. 

Proof. By repeated application of Th. 14.11 we get 




(4.4) 


where 

OO 

F^^\u:) = JJ 

k=m 

Since the symbols A^^’'^fiz) have the hat function symbol A^^'^{z) as funda¬ 
mental symbol, the inhnite product converges [TT]. Moreover, for m > 0 
the symbol A^^’^^^z) has no zeros on the unit circle, so that has no 

zeros, too. Thus, , m > 0, has zeros of the same order as 
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does, and the claim follows. 

For m = 0, 

^(i,o)(g-*f) g simple zero for oj = 2{2a + l)7r, a G Z, while pgg 

a zero of order n — 1 for a; = davr, a E Z\{0}. Thus, has a simple 

zero for oj = 2a7i, a E Z\{0}. □ 

The Strang-Fix conditions allow us to conclude that polynomials of suit¬ 
able degree are contained in the space generated by the 2“"^-translates of 
(^(n,m), space 

= spm - 2-^a) : a eZ] , m E Z+. (4.5) 

Theorem 4.3. Let 11^ be the space of polynomials up to degree d. For any 
n > 2 and m > 0, the space contains Tln- 2 , he. for any polynomial 

p E n „_2 there exists a sequence of real numbers : a G Z} such that 

p = J2 7a- 2"™a). (4.6) 

aGZ 

The space contains the space of the constants IIo. 

Note. The roots of the nonstationary symbols (I3.12p cannot fulhll the hy¬ 
potheses of Th. 1 in |32]. Thus, the rehnable functions cannot generate 

exponential polynomials. Actually, we are interest in the construction of ef- 
hcient decomposition and reconstruction formulas for general (possibly non 
exponential) signals. To this end high algebraic polynomial generation is 
more effective since it induces a high number of vanishing moments in the 
analyzing wavelet. 

From the results above it follows the approximation order of the system 

Corollary 4.4. For any n > 2 and m > 0, the system has approx¬ 

imation order n — 1, i.e. for any f E T 2 (M) there exists a constant Cf, 
independent from m, such that 

inf J|/-/^|| 2 <C; 2 -™("-i). ^ 47 ) 

The system has approximation order 1. 
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Note. Even if the system reproduces just the constants, polynomials of 

higher degree can be represented by suitable integer translates of ■). 

In fact, G and its Fourier transform has a zero of order 

n — 1 for a; = 2Q!7r, a G Z\{0}. Thus, polynomials of degree n — 2 can 
be represented by the integer shifts of the dilate and the system 

{^(n,0)(2-l . 

— a) : a G Z} has approximation order n — 1 (cf. |H]). 

Interestingly enough, the convolution property fl4.3p allows us to prove a 
differentiation rule for the functions 

Theorem 4.5. Let he the backward finite difference operator defined re¬ 
cursively as 

Vk/=i(/-/('-A)), Vj/ = /, = r>l. 

For n >2, the derivatives of are given by 

p)r ^(n,0) ^ ^ r <n-l, 

(4.8) 

^ r<n-l, m>0. 

Proof. We will prove fl4.8p just for m > 0; the case m = 0 can be proved in 
a similar way. 

From (14.21) and the hrst of (14.4p it follows 

(_D (p("-.™))(ci;) = ioj ^("■’™)(c(;) = iuj = 



The inverse Fourier transform gives 

which is the derivation rule for r = 1. Repeated applications of this rule give 
the derivation rules of higher order. □ 
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From fl4.8p and the variation diminishing property, we can infer that the 
derivative m G , has the same behavior of the hnite difference 

= meZ+. (4,9) 

Q! = 0 ' 

with h = 2~^ for m = 0 and h = 2“™ for m > 0. In particular, the number 
of strict sign changes of is not greater than the number of strict 

sign changes of the sequence 

= {c^ : 0 < a < r} = {(—!)“ : 0 < a < r} , 

so that we can infer the shape of from S~ (c^). 

Theorem 4.6. For any n > 2 and m E is bell-shaped, i.e. 

centrally symmetric, strictly increasing on [0, \supp (^^"’"^^1/2], and 
its second derivative has just 2 sign changes. 

Proof. Since any mask is centrally symmetric, any is centrally 

symmetric, too. As a consequence D is centrally antisymmetric and 

is centrally symmetric. For m > 0 fl4.8p and fl4.9p give 

D^{n,m) ^ 1 _ h)) , 

(4.10) 

JJ2 ^ _ 2 - k) + - 2 k)) , 

where h = 2“"*, while for m = 0 fl4.8p - fl4.9p together with fl3.4p yield 
D (^("’0) = 1(1 (^(^-14) (x) - 1 (^(^-14) {x - 2 h)) , 

]J2 ^{n,0) ^ ^(1 _ 1 (4.11) 

1 (^x — 2 h) + ^ (^(«^-i4) (a; — 3 h)) , 

where h = 2~^. Since is positive in (0, |supp Lptr>—i,m)y^ ysxi- 

ishes for x < 0, and are positive for 0 < x < h. Now, 

by the variation diminishing property from the hrst relations in (Id.lOp and 
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fl4.1ip it follows that has just one sign change which has to be in 

|supp thus, > 0 in ( 0 , |supp and < 0 in 

(|supp (y 9 ("'’™)|/ 2 , |supp As for the second derivative, a direct com¬ 
putation shows that jg negative in |supp so that 

has at least two sign changes; but from the variation diminishing property 
and the second relations in fId.lOp and fid.lip it follows that cannot 

have more than two sign changes, so proving the claim. □ 


5 Nonstationary Prewavelets 


From the results in the previous sections it follows that, for any n >2 held 
hx, the spaces m G Zi+, generate a nonstationary multiresolution 

analysis as dehned in Section 2. 

We note that the L 2 (M)-stabihty of the basis implies that the 

symbol of the autocorrelation of i.e. the polynomial 


(n,m) 

p)p 


(w) 


^ + 2-^a) 

aez 


itA; 2 a 


2m + , 

aeZ 


(5.1) 


is non vanishing for any a; G M. The vector = [ 77 ^"’™'^]^, where = 

Jjg (p("’™')(--l- 2 “™'Q()), is the eigenvector corresponding to the eigenvalue 

1 of the transition operator 


(y(n,m) A)„ = 2 ^ A;3 , « e Z, A G ioiZ ), (5.2) 

/3ez 

where , a G Z}, is the autocor¬ 
relation of the mask . The sequence { 77 !”’™'^} is positive, compactly 

supported with 


(n,m) 


suppir]^’”"^} 




m = 0, 


[-77 - 1, 77 + 1] , 


777 > 0 , 
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(5.3) 


and centrally symmetric. As a conseqnence g M with 

0 < < pJT’^Ho) = 2”^. 

Since 

^<,"•”■>(2“^) = E„ez (-1)“ (/r + Z-^o))) = 

(5.4) 

= Eaez * (/.("’-))) + 2—«), 

at any level m the basis is a non orthogonal basis. 

The nonstationary mnltiresolntion analysis {!/("■’”*) ; m G Z+} allows ns 
to dehne a wavelet space seqnence ; m G Z+}, where each space 

\Y(n,m) jg orthogonal complement of in Since is 

non orthogonal, orthogonal wavelets with compact snpport do not exists. On 
the other hand, dne to the L 2 (M)-stability and the compact snpport of each 
i_p(n,m) ^ it is always possible to constrnct compactly snpported semiorthogonal 
wavelets [El Th. 3.12]. 

The explicit expression of the prewavelet of minimal snpport can be ob¬ 
tained generalizing to the nonstationary case the results in |27] (see also 

HE]). 

Theorem 5.1. For any n>2, the functions 

n+1 

a=—n 

(6.6) 

n+1 

^(n,m) _ (-1)^ a) , m > 0 , 

a=—2n 

where 

gin,m) ^ f ^(n,m+l)^. ^ a) , « G Z , (5.6) 

Jr 

are the prewavelets generating the wavelet spaces 

tyO’™) = -2"™a) : a G Z} , m G Z+. (5.7) 

The functions compactly supported. Moreover, any system 

^(n,m) _ « G Z} (5.8) 

is L 2 (M.)-stable and linearly independent. Finally, vanishing 

moments. 
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Proof. A straightforward computation yields 

/ - 2"™/3) = 

Jr 

= V (-1)“ - 2-™/3) = 

a 

E / ^\a (n,m) (n,m) _ „ 

I -‘-J Sa-l 92/3-0 ~ ) 

a 

for any /? G Z, i.e. is orthogonal to the space Due to the 

support of it follows 


(T, 




[—n — 1, n] , m = 0 , 


[—2n — 1, n] , m > 0 , 
so that is compactly supported with 

[-|,n + l] , m = 0, 


(5.9) 


SUpp = 

[-2-"^u,2-’"(n +1)], m>0. 

Now, the wavelet equations fl5.5p in the Fourier space reads 

^in,m) ^ ^(n,m) (g-i^ ) ^(n,m+l) ^ m G Z+ , 

where S''^'^\z) = 2 ;", = (—1)“ An explicit calculation 

gives 

2™+V) . (5.10) 

As a consequence, 

(/« + 2“"’a)) = 

= 2”^ Eaez + 2™+%a) ^ = 

= 2”^ |A(”’™)(-e*(5^+™))|2 pSr’™+^^(a; + 2™+i7r(a + l)) 

X I (a; + 2"*+%a) f , a; G M, 


X 
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is non-vanishing. It follows that at any level m, the system gen¬ 
erating the wavelet space is linearly independent and L 2 (M)-stable. 

Moreover, is non orthogonal with 

[ + 2-^a) = ^ ^ ’ 

»/ IR. ^ 

where — {Xl, *} autocorrelation of the sequence 

{pi”’”*}. 

Finally, since 1_ from Th. I4.3l it follows that = 

0, 0 < d < n — 2, so concluding the proof. □ 


Theorem 5.2. For any n > 2 and m G is the unique minimally 

supported wavelet generating the wavelet space 

Proof. The orthogonality conditions f 'ip^'^’'^\x) — 2~'^/3)dx = 0, 

Jr 

/3 G Z, can be proved to be equivalent to the conditions 

dgi = 0, /3 gZ. (5.11) 


Let d!i'''^\z) = 4a d!jj''^\z) = Y,a 4a+i similarly, for the 

polynomial^ 2 ;", let = 


Zlo so that = d!i'’'^\z‘^) -|- z d^o''^\z‘^), and g^'^’'^\z) = 

g{ri,m)^^ 2 ) _|_ ^ g{ri,m)^^ 2 y foUows that couditious fIS.lip are equivalent to 


+zd^^'^\z) g^^'^\z) = 0 


(5.12) 


Thus, is the minimally supported prewavelet if and only if is the 

minimally supported polynomial satisfying fl5.12p . i.e. if and only if d!i^'^\z) 
and d^’'^\z) have no common zeros or, equivalently, and d‘^”’™'^(—z) 

have no common zeros for 2 ; = e* 2 ™-+i. But this easily follows from fl5.10p since 
any symbol has only zeros with negative real part and 

is positive. □ 
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We note that is the 2™-dilate of the prewavelet constructed in [T^ 

and Th. 15.11 and Th. 15.21 generalize to the nonstationary case the stationary 
prewavelet constructed in [27] . 

Even if the prewavelet coefficients do not have an explicit expres¬ 

sion, they can be evaluated efficiently by an iterative algorithm. 

Theorem 5.3. Let = [g^a''^\a G be the nonstationary pre¬ 
wavelet coefficients, and let = [g^\a G be the prewavelet coef¬ 

ficients corresponding to the stationary fundamental mask . 

For any n>2 and m G held fixed, consider the iterative procedure 


r pin,m) ^ ^ 

1 p(n,m) _ p<(n,m) p{n,m+l) 

k > 0, 

(5.13) 

where with 



(n,m) _ An,m) (n,m+l) 

“/3 “a+2/3 ’ 

p e , 

(5.14) 

and 

(n,m) j [ 2, 77, -|- 1] , 

“ \ [-2(n-l-l),n + l] , 

o o 

II A 

(5.15) 

The sequence converges strongly to 

the following error estimate holds 

when k ^ oo. 

Moreover, 


_ j^(n,m+A:)|| ^ 

(5.16) 


where 'yn,m is a positive constant independent from k. 

Proof. Using the refinable equation fl3.4p in fl5.6l) we get C 2 a-p 

which in matrix form can be written as Repeated 

applications of both algorithm fl5.13p and the relation above give 

< II nf=o (-Y(n,m+0|| . _ j\^(n,m+fc+l)|| _ 
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Now, let where = X] ®i"f 2 / 3 - piC^"^^) = 

1. Since is the fundamental mask of the mask sequence {a^”’™'^}, it 
follows EmGZ+ < oo, and 


m+fc 


Yl < exp I ^ - C'(")|| I < cx), 

m€Z_|_ 


l=m 


(cf. [m Prop. 2.1]). Moreover, limfc_,.oo , thus. 


lim = 0 , 

k^oo 

and the claim follows with = exp (j2mGZ+ 


□ 


6 Nonstationary Biorthogonal Bases 

The rehnable functions (^("■’”*1 and the prewavelets are linearly inde¬ 

pendent and L 2 (®)-stable, but they are not orthogonal. As a consequence, 
the dual bases of and in and respectively, have 

infinite support. Nevertheless, compactly supported bases giving rise to ef¬ 
ficient reconstruction and decomposition formula of a given discrete signal, 
can be obtained introducing biorthogonal bases. 

The theory of biorthogonal bases in the stationary case [3] can be gen¬ 
eralized to the nonstationary framework (cf. [22l |2H [32]). For any mul- 
tiresolution analysis ; m G Z+} we can introduce a biorthogonal 

multiresolution analysis ; 777 , ^ with 

yin,m) ^ Y{n,m+1) ^ m ^ Z+ , (6.1) 

and biorthogonal wavelet spaces : m G Z+} and {ITA’™-) : m G Z+}, 

such that for any m G Z+ 

pp(n,m) _ Y{n,m+1) q ^(n,m) ^(n,m) _ -^(n.m+l) q ^(n,m) 

( 6 - 2 ) 


At each level m, the spaces pA’™), ppl".™-) and are generated 

by the 2 “™-integer translates of the biorthogonal functions 
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^(n,m) ^ respectively, satisfying the following biorthogonality con¬ 

ditions: 




■ - 2 -^( 3 )) = 



•- 2 —/ 3 )) = 5 „^ 



•- 2 —/?)) = 0 , 



■- 2 -^ 13 )) = 0 . 


We stress that all the biorthogonal functions and 

m G cannot be obtained each other by dilation. Thus, none of the 
biorthogonal spaces at level m is a scaled versions of the spaces at level 0 . 

Biorthogonal bases for the exponential splines were constructed in |32]. 
Here, we want to construct the biorthogonal bases associated with the non¬ 
stationary rehnable functions m G Z_|_. 

As a consequence of (16.1|1 and (I 6 . 2 p . the wavelet belongs to 

while the biorthogonal functions and belong to so that 



_ 2-(”^+i)«), 

m G Z+ , 



(• - 2-(™+i)a), 

m G , 

(6.4) 



m G Z+ . 


Moreover, perfect reconstruction at each level m is 
m G the biorthogonal symbols 

guaranteed 

if for any 



(6.5) 

satisfy 

-z ^) = 1 , 

( 6 . 6 ) 

with 







(6.7) 


Identity (16.611 is a Bezout’s equation which has a unique polynomial solution 
j^in,m) given degree [ 7 ]. 
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Under mild conditions on the symbols and biorthogonality 

conditions fl6.3p guarantee that the biorthogonal bases 

^ _ 2-™a), tt G Z} , - 2-™a), a G Z} , 

^ {(^(n,m) j. _ 2-™a), tt G Z} , - 2-™a), a G Z} , 

( 6 . 8 ) 

are L 2 (M)-stable, so that, for any / G L 2 (M), the following decomposition 
formula holds 


/ = /»..+ E E (/• - 2-”‘a), (6.9) 

m>mo aGZ 


where 

/mo = ( 6 . 10 ) 

is the mo-level approximation. Some examples of nonstationary biorthogonal 
bases will be given in the next section. 

The sequences {oq,}, {oq}, {qa}, {Qo} are pairs of biorthogonal FIR hlters 
that give rise to the decomposition and reconstruction algorithms 


\ m 


—y 

\pi ^ 


~{n,m) Am+1 
“^-2a ^/3 5 


/-m _ 1 ~U,m) ,m+l 

8a “ AT ^/9-2a ^/3 > 

/3ez 


( 6 . 11 ) 


^m +1 


1 

7! 


E ^(»i,m) I „U.m) 

®q:- 2/3 ^/3 + '?a-2/3 8/3 

. /3gZ /3gZ 


( 6 . 12 ) 


which can be efficiently used for the analysis and synthesis of a given data 
sequence A° = {A° : a G Z}. 


7 A case study 

In this section we give some examples of both nonstationary prewavelets and 
biorthogonal bases in the case when n = 3. In this case the nonstationary 
refinable functions m > 0, belong to C'^(R), i.e. they have the same 

smoothness as the cubic B-spline. Interestingly enough, any with 

m > 0 has the same support as j_g_ j'g^q . 2“"*], while is more 

localized in the scale-time plane having suppcp^^’^^ = [0, 5/2], a property that 
appears very useful in applications (see the example below). In order to 
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Table 1: Numerical values (rounded to the forth digit) of the mask coefficients 


(3,m) (3,m) 

% , b 

^(3,m) ^ 

...,8. 

Here fi = 

1.1 




m 

0 

1 

2 

3 

4 

5 

6 

7 

8 

(3,m) 

% 

0.5 

0.0313 

0.0452 

0.0508 

0.0537 

0.0555 

0.0567 

0.0576 

0.0583 

(3,m) 

0 ) 

0.5 

0.2500 

0.2500 

0.2500 

0.2500 

0.2500 

0.2500 

0.2500 

0.2500 

(3,m) 

^2 

0 

0.4375 

0.4095 

0.3984 

0.3925 

0.3889 

0.3865 

0.3848 

0.3835 




Figure 1: The nonstationary mask coefficients listed in Tab. [U (left) and 
(right). The stationary mask of the cubic B-spline and the cubic B-spline 
itself are also displayed (dashed line) 


obtain refinable beses and nonstationary filters significantly different from 
those ones generated by the cubic B-spline, we choose /i = 1.1 as a value for 


the tension parameter. 

The coefficients of the mask = 
while the coefficients of the mask 


(3,0) 
0 i 


{a, 

= {fli 




(3,m) 
0 i 


(3,m) 

a\ 


(3,0) 

are % 

(3,m) 
^2 • 


= 


(3,m) 


1 

2 ’ 


(3,m)^ 

“4 


for m > 0 are 


(3,m) _ (3,m) _ (3,m) _ (3,m) _ J- (3,m) _ J- _ 

Uq — (44 — Z, , U4 — Us ~ ~ 2 ^ 


The numerical values (rounded to the forth digit) of the mask coefficients 
are listed in Tab. [H while their behavior is shown in Fig. [H (right). The 
behavior of in comparison with is displayed in Fig. [1] (left). The 
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Figure 2: The prewavelet mask { 5 '*'^’°^} (left) and (right) 


nonstationary prewavelets are given by: 

V.(3,»)= Y. (-ir9h"T<“’'>(--2-'a). 

o=-3 

4 

^(3,m) _ (-I)“^2l? - 2-(”^+^)a), m > 0 , 

a =—6 

where the prewavelet coefficients {ga’"'^} can be evaluated by the algorithm 
in Th. 15.31 From fl5.9p it follows that = [—3/2,4], while for m > 0 

= [—2“™'3, 2“™ 4], We notice that is more localized in 

the scale-time plane than both ■^(3’™-), m > 0, and the B-spline prewavelet. 
The prewavelet mask coefficients, rounded to the forth digit, are = 

= -0.0015, g^^i^ = = 0.0259, g^^^^ = -gf’^^ = -0.1479, 

£3.0) ^ _^(3,o) ^ 0.3244. 

In Fig. [2]the behavior of and .0(^.0) are displayed. Finally, we give the 

explicit expression of the biorthogonal masks m > 0. It is well known 

that the biorthogonal mask of a^^.o) is 5 (^, 0 ) _ /£3,o)^~(3,o)| _ 

In order to fulhll conditions ensuring the existence of the biorthogonal re- 
hnable function £3.™-), for m > 0 we construct the biorthogonal mask a^^.™') 
with support [0,14]. Its explicit expression is given by 
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(3,m) _~(3,m) 


Og = a 


14 


~(3,m) ~(3,m) 

' = a ' - 


13 


(3,m) _~(3,m) 


Clr\ - Ct 


12 


(3,m) _~(3,m) 


do — d 


11 


(3,m) _~(3,m) 


a;, = a 


10 


~{3,m) ~(3,m) 

= On 


~{3,m) ~(3,m) 

% = «8 


07 = 


^ —4—h 


-4 + 2 ^^ 


= -j-(128 + 2^+’^ + 5 • 4^+^ + 5-8'") 

—4 + 2^ 

= — - -j-(128 + 2®+^ + 5 • 4}+’^ + 5-8'^) 

—4 + 2" 

Q—3—h 

= -^-^(640 + 7 • 2®+'* + 33 • 4^+^* + 29 • 8^ - 5 • 

—4 + 2” 

2-9-2-h 

= -^(128 + 3 • 2^+^* + 17 • 4^+^" + 17 • 8'*) 

-4 + 2^^ ' 

o—3—h 

= — -7-(1152 + 15 • 2®+'* + 133 • 4}+^ + 89 • S'* - 39 • IG'*) 

—4 + 2^ 

= — -^(128 + 2®+'* - 123 • 4^+^* - 123 • S'*) 

-4 + 2^^ ' 

Q—3—h 

= -^(G40 + 9 • 2®+'* - 81 • 2^+*^-'* + 105 • 4^+'* + 577 • S'*) 

-4 + 2^* ^ ' 

(128 + 3 • 2®+'* + 81 • 4^+'* - 175 • S'*) 


where h = 3 + m The biorthogonal mask coefficients (rounded to 

the forth digit) are listed in Tab. |2l In Fig. [3] the behavior of and 

~(3,o) displayed. The biorthogonal wavelet mask coefficients q(^’”*) and 
q(3,rr*) obtained by = (—l)"a^a™\, = —(—l)“a?a™\. In 

Fig. mthe behavior of and is displayed. 

Just to show how the properties of the constructed nonstationary biorthog¬ 
onal filters can affect the analysis of a given signal, we evaluate the coefficients 
{A™} and {C^}, obtained after three steps of the decomposition algorithm 
flb.lip . when the starting sequence is a spike-like signal (see Fig. [5] (left)). 
The coefficients are plotted in Fig. [5] (right) in comparison with the coeffi¬ 
cients obtained when using the stationary cubic spline biorthogonal hlters. 
The hgure shows that the nonstationary decomposition algorithm has higher 
compression properties: actually the number of nonzero coefficients are 26 in 
the nonstationary case, while they are 39 in the B-spline case. 
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Table 2: Numerical values (rounded to the forth digit) of the mask coefficients 


~{3,m) r 

Qa , tor a 

= 0 ,... 

, 7, and m 

= 0 ,.. 

.,8 





m 

0 

1 

2 

3 

4 

5 

6 

7 

8 

«o 

0.5 

0.0011 

0.0021 

0.0026 

0.0030 

0.0064 

0.0034 

0.0035 

0.0036 

•~(3,m) 

«i 

0.5 

- 0.0085 

- 0.0114 

- 0.0129 

- 0.0138 

- 0.0288 

- 0.0148 

- 0.0151 

- 0.0154 

•~(3,m) 

0 

0.0066 

0.0028 

0.0005 

- 0.0010 

- 0.0039 

- 0.0027 

- 0.0032 

- 0.0036 

~(3,m) 

^3 

0 

0.0574 

0.0760 

0.0857 

0.0914 

0.1905 

0.0979 

0.0999 

0.1014 

~(3,m) 

U4 

0 

- 0.0810 

- 0.0790 

- 0.0768 

- 0.0752 

- 0.1480 

- 0.0732 

- 0.0725 

- 0.0720 

•~(3,m) 

«5 

0 

- 0.1998 

- 0.5108 

- 0.2834 

- 0.2998 

- 0.6211 

- 0.3180 

- 0.3236 

- 0.3278 

•~(3,m) 

0 

0.3233 

0.3241 

0.3237 

0.3232 

0.6456 

0.3225 

0.3222 

0.3220 

(JjiJ 

0 

0.8019 

0.8816 

0.9212 

0.9443 

1.9187 

0.9698 

0.9776 

0.9835 




Figure 3: The hrst 8 nonstationary biorthogonal masks (left) and 

(right) 




Figure 4: Graphs of (left) and -0(3’°) (right) 
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0.9 - 
0.8 - 
0.7 • 

0.6 • 

0.5 • 

0.4 - 
0.3 - 
0.2 - 
0.1 • 

fll-'-■-'-U-■-'-■- 

0 200 400 600 800 1 000 1200 1 400 1600 



Figure 5: The spike-like initial sequence (left) and the decomposition coef¬ 
ficients obtained after 3 steps of the nonstationary decomposition algorithm 
(upper right). The decomposition coefficients obtained by the stationary 
cubic biorthogonal filters are also shown (bottom right) 

8 Conclusion 

We studied the properties of a class of refinable ripplets associated with se¬ 
quences of nonstationary scaling masks. One of the most interesting property 
of these functions is in that they have a smaller support than the stationary 
rehnable ripplets with the same smoothness. This localization property is 
crucial in several applications, from geometric modeling to signal processing. 
After proving some approximation properties, such as Strang-Fix conditions, 
polynomial reproduction and approximation order, we proved also that any 
rehnable function in the family is bell-shaped, so that they can efficiently 
approximate a Gaussian. Moreover, since these rehnable functions gener¬ 
ate nonstationary multiresolution analyses, we constructed the minimally 
supported nonstationary prewavelets and proved that their 2^-shifts form 
L 2 -stable bases. We note that this construction can be generalized to other 
class of nonstationary rehnable functions, like exponential splines. Moreover, 
we constructed nonstationary biorthogonal bases and hlters to be used in ef- 
hcient decomposition and reconstruction algorithms. 

The localization property of the rehnable ripplets we studied implies that the 
corresponding nonstationary wavelets have a small support too, a property 
which is very desirable in the case when the relevant information of a func¬ 
tion to be approximated or of a signal to be analyzed are focused in small 
regions of the scale-time plane. The preliminary test in Section 6 shows 
the good performances of the constructed nonstationary wavelets in a simple 
compression test. More tests will be the subject of a forthcoming paper. 
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